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Based on an approach introduced by Gerjuoy, Rau, and Spruch, we con- 
struct variational principles in a systematic way for the nonlinear Schrodinger 
equation and obtain new variational principles for the case of Ginzburg- 
Pitaevskii-Gross equation (GPG) which is believed to describe accurately 
the Bose-Einstein condensation at zero temperature. As an application of 
these variational methods, a variational iteration method is proposed for 
calculating eigenvalue (chemical potential) and wave function for the GPG 
equation 
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1. Introduction 

The concept of the Bose-Einstein condensation (BEC) [1] has been known 
for 73 years, and has been used to describe all physical scales, including liquid 
A He, excitons in semiconductors, pions and kaons in dense nuclear matter 
(neutron stars, supernovae), and elementary particles [2]. It is only a few 
years ago that the BEC phenomenon was observed directly in dilute vapors 
of alkali atoms, such as rubidium [3], lithium [4], and sodium [5], confined in 
magnetic traps and cooled down to nanokelvin temperatures. 

These remarkable experimental observations [3-5] have stimulated much 
activities in experiments [6-11] and in the theoretical investigation of inho- 
mogenous Bose gases. It is impractical to quote all of numerous theoretical 
papers published during the past few years. We quote here only a frac- 
tion [12-28] of the theoretical papers, which are relevant to this paper and 
are based mostly on the Ginsburg-Pitaevskii-Gross (GPG) equation (also 
known as Pitaevskii-Gross (PG) equation) [29-32], which is believed to de- 
scribe accurately the experimental results well below the onset of BEC. The 
GPG equation is a nonlinear Schrodinger equation (NLSE), which is formally 
identical to a Ginzburg-Landau type of field equations (the Ginzburg-Landau 
theory is also known as ip theory [36]). The NLSE is one of the nonlinear 
differential equations, the study of which led to fundamental advances in 
nonlinear dynamics [37]. 



Since the vapors used in the experiment are dilute, the average distance 
between the atoms is substantially larger than the range of the interaction. 
Hence, the Bose gas is expected to be dominated by two-body collisions, 
which can be well described by the s-wave scattering length. Thus, the 
Ginsburg-Pitaevskii-Gross theory developed for weakly interacting bosons is 
ideally suited to describe BEC observed in the experiments. 

The numerical methods based on variational principles (VP) are the most 
effective approximation method and have been used to solve the GPG equa- 
tion using the numerical minimization of the energy functional [14, 16, 21, 
25, 28]. However, it is expected that the wave function determined by the 
variational method is much less accurate than the stationary value of the 
functional itself (e.g. energy), because the later value is only indirectly re- 
lated to the wave function and depends weakly on the trial wave function 
near the stationary point. We note in contrast that, for the case of the linear 
Schrodinger equation, one can obtain the criteria of accuracy not only for sta- 
tionary quantities but also for the wave function using Schwinger variational 
principles [35]. 

In this paper, we present a systematic construction of VP which is free of 
the above shortcomings, based on the approach proposed by Gerjuoy, Rau, 
and Spruch (GRS) [36]. We introduce a functional whose optimal value is 
the wave function at given values of its arguments. In section 2, we give a 
brief review of the GPG (or PG) equation as it relates to BEC. In section 
3, we give an outline of GRS approach for variational principles. Section 
4 describes our derivation of a functional, stationary value of which is the 
chemical potential, and section 5 describes the VP for the wave function. In 
section 6, we describe our formulation of iteration method based on our VP 
for the wave function. Summary and conclusions are given in section 7. In 
the Appendix we describe an example of the application of our variational 
iteration method to a nonlinear integral equation. 

2. Nonlinear Schrodinger Equation 

For a dilute gas of bosons at zero temperature mean field theory yields 
a nonlinear Schrodinger equation (NLSE) for the condensate wave function 
[32]. The condensate wave function (or order parameter) may be described 
by a self-consistent mean field NLSE known as the Ginsburg-Pitaevskii-Gross 
(GPG) equation (also known as Gross-Pitaevskii (GP) equation [29-31]. For 
condensed neutral atoms in an external trap, the GPG equation takes a form 
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VV(r, t) + V trap (r)^(f, t) + V H (f)7P(r, t) 
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where ip(r,t) is the BEC wavefunction,m is the atom mass, and Vt rap (r) 
is a trap potential. Vff(f) is the mean Hartree interaction energy for the 
condensate 

V H (r) = g N \il>(r,t)\ 2 , (2) 
where N is the total m number of the BEC atoms, and 

<iiih 2 a . . 

9o = 3 

m 

with the scattering length a. To obtain a stationary solution, we insert 
ip(r,t) = exp(—i l f-t)^(r) (where /x is the chemical potential) into Eq. (1) 
to find the following equation for ip{f) 

h 2 

-^V 2 V(r) + V trap (f)^(f) + g N \^(f)\ 2 ^(f) = /i ^(f), (4) 



df\^j(f)\ 2 = l. (5) 

For an anisotropic harmonic oscillator trap, V trap (f) is given by 

V trap (f) = m(uj 2 x x 2 + uj 2 y 2 + oj 2 z 2 )/2. (6) 

In this paper, we develop the variational principles, for solving Eq. (4) based 
on the approach proposed by Gerjuoy, Rau, and Spruch [36]. The GRS 
variational principles for the time-dependent Eq. (1) will be given elsewhere. 



3. Variational Principles of Gerjuoy, Rau, and Spruch. 

In this section, we give a brief outline of the unified formulation of the 
construction of variational principles, developed by Gerjuoy, Rau, and Spruch 
GRS [36]. 

We start with the following set of equations 

B i (<j>) = 0, 2 = 1,2,..., (7) 

where designates a set of function are precisely defined by relations (7). 

Now, let 0*) be a trial estimate of exact but unknown (0, 0*), and let 
(L\, L\), be a trial estimate of a (L\ U) to be specified. Then, the functional 
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[F] = F v ((f> t , 4>*i L, L) is a variational functional for the desired F((f>, 0*) which 
we seek, if 

F v (<f>,<j>*,L t ,L t )=F( ( f ) , ( f ) *), (8) 

and 

S[F] = F v {4> t , L t , L t ) - F(4>, 0*) = O(o0 2 , 5L\ o0oZ), (9) 

where 

F = Ll + 5L\ U = U + 5L\ = t + 8<j>, 0* = <f>* t + 5(f)*. (10) 
Now, consider 

[F] = F(0 t , t *) + Y,m\Bi{<t> t )) + (B^Ll)). (11) 

i 

Using Eqs. (9 - 11) we can determine the necessary relations specifying 
(L*,L l ), whose estimates (F,F) appear in Eq. (9). The essential point is 
that there are no restrictions on the nature of 0,F(0), and -8(0). There- 
fore the GRS method is a general method for construction of VP for linear, 
nonlinear, differential, and integral equations. 

4. Variational Principle for the Chemical Potential 

In this section we derive a functional whose stationary value is the chem- 
ical potential /x - We start with Eq. (4) 

(T + V trap + g N \^\ 2 )tl) = (4a) 

where T is the kinetic energy operator 

h 2 

T = -^V 2 . (12) 
2m y ' 



Now, let us consider a functinal [«o] 

N = ii t + \%{M^t)-l) + {Lt\(T + V trap + g N^ t \ 2 -^ t )\^ 
+ ({T + Vtra V + 9oN \A\ 2 ~ lk)i/>t\L t )], 

where 



(13) 



A * (lm + d) t \L t y (14) 
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and L t (r) is solution of linear equation 

(T + V trap + 2g N \^ t \ 2 )L t + g N $L* t - (i t L t = (15) 
with boundary conditions 

lim L t (r) -> 0. (16) 



r— >oo 



To prove that [/to], Eq. (13) is a stationary expression, let us calculate 5[/io\, 
let 

ipt = -0 + <ty>, ^* = ^* + <ty*, 
^ = A*o + 8fx, L t = L + 5L. 



(17) 



We also assume that 

Sfj. = 0(6if>), 5L = 0(#). (18) 
Substitution of Eq. (17) into Eq. (13) gives 

%o] = W-A(W) + ML>)]+ 
A((<ty#> + (Vl#) + (^l(T + ^ rap + 2^ Ar |V| 2 -/i o )|L) + 

(<^% AW> 2 |L*> + ^|(T + \4 mp + 2f/ iVo^| 2 -/io)^|^) + 
{g N ^L*\6r/;)) + 0(# 2 ). 

(19) 

Using 

A = (<LhA> + <^|L>)-\ 

and Eq. (18), we obtain 

% ]=W 2 )- 
Since homogeneous equation 

(T + Kra P + 2g N \^\ 2 )L + (? iVo^ 2 ^ - \iL = (21) 

has a solution 

L = ic^, (22) 

where c is a real constant, and there is no homogeneous solution in general 
of 

(T + V trap + 2g N \^ t \ 2 )L t + g N ^ 2 L* t - &L t = 0, (23) 

our assumption given by Eq. (18), 5L = 0(5ip), has to be justified. In 
general, we need to modify Eq. (15). (For the case of the linear Schrodinger 
equation, see [36]). However, for ground state of the GPG equation, we do 
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not need to modify Eq. (15). In this case, ipl = ip t , and from Eq. (15) we 
write 

(T + V trap + 2g N \^ t \ 2 )R e L t + g N \^ t \ 2 R e L t - fi t R e L t = -iP t , , . 
(T + V trap + 2g N \ij t \ 2 )ImL t -g N \ij t \ 2 ImL t -fi t ImL t = 0, ^ 

where 

L t = ReL t + HmL t . 

For this case, we need to solve Eq. (15) for real solution L t (L* t = L t ) and 
hence the problem will not arise. 

5. Variational principle for wave function 

In this section, we attempt to solve the following problem how to find 
variational estimates of the BEC wave function. To do this, we introduce a 
functionals [Reip]{f) and [Imip](r) whose stationar values are [Reip]{r) and 
Imipif) respectively 

[Re^](f) = l^ t ^+^^)-ljLl\r,n[T+V trap + g N \ij t \ 2 

- ^(^dr" - ± f[(T + V trap + g N \^ t \ 2 - ^ t )Mn}*LW^dr>, 

(25) 



[Im^]{r) = UM?) ~ - ^[/Lf (r,F)[T + V trap + <7oWt| 2 - IM^lf)*? 

+ I J[(T + V trap + U N \ij t \ 2 - ^)^(r^)]*LKr, f-')df], 

(26) 

where L\(r, f), and L R (f, r") are solutions of linear equations 

[T + V trap + 2g N \ij t \ 2 - f i t ]L t R (ry)+g N ij 2 L t ;(ry) =S(f-f), (27) 



[T + V trap + 2g N \^ t \ 2 - nAL* (f, 7) - g N ^ 2 Lf (f, 7) = S(r- f). (28) 
Varying Eqs. (25) and (26) we obtain 

5[Re^\{r) = J dr , ^{6(f-f)(6i/)(f)+Si/)*(f))- 

\[{[T + V trap + 2^ V |^| 2 - ^ t \L R (f,r'))*5iP(r > ) 
+ g N ^* 2 (f)L R (f,f)Sip(f) (29) 
+ [T + V trap + 2g N \^ t \ 2 - Vt]L R (r, 7*)<ty V) 

+ goNonnL^^rin]} + w 2 ), 
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and 

5[Irr4\{?) = Jdr*{±{5(r-r*)(5ij(r*)-5ij*(r*)) 

- Ud T + v trap + 2<?o^cM| 2 - Ai t ]LKr,rW(r) 

- g N ij* 2 (7)^(^7)5^) (30) 

- [T + V + 2g N \^\ 2 - ^L^^rin 

+ g Q N Q ^\f')L 1 f{r,f')Srj}*{T f )\\ + 0(5ip 2 ), 



where we have assumed 

fJ>t = {J>o + Sfx,Sfx = 0(5ip 2 ), (31) 

and 

L\, = L RJ + SL RJ , 5L rj = 0(SjP). (32) 
Substitution of Eqs. (27), (28) into Eqs. (29) and (30) leads to equations 

5[Re4>](r) = 0(5^), mi 

8[Im4(r) = 0(# 2 ), 1 J 

which prove that Eqs. (25) and (26) are variational functionals. 

Using Eqs. (25) and (26) we can derive a variational principle for if) itself 

M (r ) = [Re^j] (f ) + i [Im^j] (f) = ^ (?) + / R% (r, 7) [T + V trap 

+ 9oN \ip t \ 2 ~ fk]ipt{r')dr' (34) 
+ / l(T + V trap + g N \^ t \ 2 -fi t )MnrR\(r,ndr>, 

where 

R\{r,7) = -\(L R (ry) - L\(r,r% (35) 

R\{ry) = -\{L R {ry) + L\{r,r')). (36) 
It is convenient to rewrite Eq. (34) as follows 

M(0 = -2 / ^(0^o|^(OI 2 ^(r,0^ ^ 

- 2 / R%(ry)g N \4> t (i*)\ 2 M?)dr*. [6<) 

Three comments are appropriate here. First, for the ground state of the GPG 
equation, we have ip* = ipt, and from Eq. (21) 

= fMn9oN Mn\ 2 L R yy)dr> , , 

+ fL»(ry)g N \'il>t(r')\ 2 Mr')dr'- 1 J 
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Secondly, if Ln(f,r') = L'j t (f,r f ), the following equation 

[T + V trap + 2g N \^j\ 2 - fi ]L R + g N ^ 2 L R = 5(f- f) (39) 

has a unique solution, since in this case the homogeneous equation has only 
a trivial solution. 

Thirdly, if ip* ^ ipt (vortex states), Eqs. (27) and (28) are meaningless. To 
see this, we use Eqs. (27), (28), (35), and (36) to obtain 

(T + V trap + 2g N \^\ 2 - f i )R 2 + g N ^ 2 R* 1 = -5(f-r'), , , 
(T + V trap + 2g N \^\ 2 - f , )R 1 + g N ^ 2 R* 2 = 0, ^ 

where ip is exact solution of the GPG equation (4). Multiplying Eqs. (40) 
on the left by ip*^) and integrating, we find, using Eq. (4), that 

/ goNomnwnnMr, n + ^rk?, o)^ = -i>(f), ^ ^ 

/ g N \ij(f')\ 2 (ij*(/)R 1 (f,/) +^(/)R*(f, r>))dr> = 0. { ^ L) 

Eqs. (41) do not have solutions if ip(r) ^ 0. As noted in [36], when construct- 
ing variational principles for certain quantities (as, for example, for ip(r)), it 
is expected that some specification of the phase of the wave function must 
be imposed; otherwise, various ambiguities and contradictions can occur (see 
Eqs. (40)). There are various ways of specifiying the phase of ip. The sim- 
plest and useful procedure is to fix the phase of ip relative to some arbitrary 
known function x( r ) through a restriction of the sort that < ip\x > is either 
purely real or purely imaginary [36]. 

Therefore we write 

bP}(r) = Mr) + lt(r)(< ip t \x >-< x\i>t >) + 

+ / Ri*(f, f) [T + V trap + g N \M 2 ~ Vt]A(7)df*+ (42) 
+ J[(T + V trap + goN \ij t \ 2 ~ Ht)Mn]*R\(r, ?)df, 



where 



(T + Vtrap + 2g N \ip t \ 2 -^t)R t 2 *(f,f)+g N ^ 2 R[(r,r') + 



(43) 



(T + V trap + 2g N \M 2 ' ^t)R{(r,n + 9oN iJ 2 Rf,(f,n + 
+ lt(r)x(f) = °> 



(44) 

7t(r) = 0.5^(f)/ < Vt|x >, (45) 
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and 

< ip t \x >=< xWt > • (46) 
Varying Eq. (42) we obtain 

^](f) = 0(^ 2 ), (47) 

where we have assumed Eq. (31) and 

SRi = 0(5ip),i = 1,2 (48) 

Using Eqs. (43)- (46) we can rewrite Eq. (42) as Eq. (37). 
6. Variational-iteration method 

To solve nonlinear problems, one frequently uses the method of iteration, 
several iteration schemes are possible [15, 37]. 

Now we formulate iteration method which directly follows from varia- 
tional methods given by Eqs. (37)and (38), 



^W(f) = / df'Q[ N - 1 \f 1 r')^ N - l \r , )df' 

where, for example, for ground state (ip( N ^ = ipW*), 

Q[ N - l \r,r>) = Qf _1) (rX) 

= goNo^-VWLg-Vftr*). 

L^~ l \r,r') is a solution of the linear equation 



(49) 



(50) 



[T+v trap + 2^i^^r-^]4^" 1) +^o(^ (jv - i) ) 2 4" r " 1) 

(51) 



where 



^ = + \(N-Q(LW-V\T + V trap + g N \^-^ - l4 If ~ 1) \^ N -% 

A^- 1 ) = l/(L^\^ N -V), 



(52) 



with 



= (^IT + V^ + V^l^), (53) 
= goNo^-^W (54) 
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and L,( N ^ is a quadratically integrable solution of the linear nonhomoge- 
neous equation 

(T + V trap + 3Vf - = -^ {N - 1] - (55) 

Eqs. (50-55) can be easily generalized to include vortex states. We note 
that one of the primary motivations of the GPG theory was the study of 
vortex states in weakly interacting bosons [30, 31]. We expect a super-fast 
convergence of the iteration process describing by Eqs. (49-55). It follows 
from the fact that, if (^ - V (JV_1) ) = 0(#), then (if, - V> (JV) )) = 0(5^ 2 ). 

To simplify the iteration process given by Eqs. (49-55) let us introduce 

^W(f) = / dr'Qf\r,r')^ N - l \r')dr' , , 

+ SdFQf \r,r'){^ N - 1 \r'))*dr'. 1 J 

This scheme is much simpler than Eqs. (49 - 55) since we need to solve a set 
of linear equations (51) and (55) only once (for N = 1). If function ip^ is 

chosen sufficiently close to ip, the Q\ N ^(r, r') and Qf'\r,f') will differ only 
by small amount. This provides the basis for the iteration process given by 
Eq. (56) [36]. To solve a set of linear equations, Eqs. (51) and (55), we can 
expand solutions in a basis of anisotropic trap eigenfunctions [15]. 

L ( t 1] (r,n = EiS^7g"" 1 V* 4 (a:)0 w (y)^W^(xO^(y / )^(^), 

(57) 

where <f> ni (y) is a one-dimensional harmonic-oscillator eigenfunction. The 
coefficients, atf 1 ' 1 ^, (3\ N ~ l \ and ^if~ l \ can be expressed in terms of 

/+oo 
dy4>m(y)4> nj (y)4> nk (y)4> ni (y), (58) 
-oo 

which can be evaluated analytically [38]. 
7. Summary and Conclusions 

Based on the GRS approach [36], we have derived variational principles 
for the nonlinear Schrodinger equation in a systematic way. In particular, 
we have obtained variational principle for eigenvalues (chemical potential) 
and wave function of the GPG equation. To the best of our knowledge, 
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these variational principles for the GPG equation are new. Using these vari- 
ational principles, we have formulated a variational iteration method, which 
is expected to have a very fast convergence. 

We are planning to generalize our results to the case of nonstationary 
GPG equation and to carry out relevant numerical calculations based on our 
method. 
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Appendix 

Our variational iteration method developed in this paper is general and can 
be applied to both nonlinear differential equation and nonlinear integral equa- 
tion. In this Appendix, we describe one example of the application of our 
variational iteration method to non-linear integral equation 

rb 

y{x) = K(x,s,y(s))ds, (x,s)e[a,b], (A.l) 

J a 

where K(x, s, y) is a continuous function, with all its variables continuous, 
and has continuous low order derivatives. We introduce a functional [y](x) 
whose stationary value is y(x) 



{x)= / K(x,t,y t (t))dt+ / L(x,x')[y t (x')- / K(x' ,t,y t (t))dt]dx' , 

J a J a J a 

(A2) 

where L(x, x') is solution of a linear integral equation 

L(M) - tds'H,,sA,t,y t (t)) = '■*<'» . (A3) 

Jo, oy t oy t 

Iteration procedure which follows directly from the variational functional 
(A.2) is 

yW{x) = S b a K{x,t,y^\t))dt 

+ $ b a L^\x,s)[y^\s) {A A) 

- ! b a K{s^ N - l \t))dt}ds, 

L(N - i)(s t) _ [ b dx > L (n-i )(s xl) dK(x',t,yM) _ dK(s,t,y^(t)) 

(A.5) 

The above iteraction procedure is equivalent to the Kantorivich-Newton pro- 
cess (KN) [37] which has developed for nonlinear integral equations: 



where 



y (N \x)=y (N ~ 1 \x) + ( j )N ^(x), (A6) 

<Pn-i{x) = £n-i{x) + J dy(N-i) <f>N-i{s)ds, {A.l) 

e N ^(x) = [ b K(x,s,y( N - l \s))ds-y( N - l \x). (A8) 

J a 
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To show the equivalence between the KN method and our method for 
solving the nonlinear integral equation (A.l), we rewrite solution of Eq. (A. 7) 

as 



<j> N - 1 (x)=f r^- 1 \x,s)e N ^(s)ds, (A9) 

J a 



where 



r<*-W) = dK(x ^- 1)(s) r^-')(.,«)«b. (Aio) 

Substitution of Eq. (A. 9) into Eq. (A. 6) gives 

yW(x) = ^K(x,t,y^(t))dt 

+ fa l-ds> dKix i^ is ' ]) (s',s)} (y(»-V (s) - (All) 

- ! b a K(s,t,y( N ^(t)dt)ds 



Using Eq. (A. 10), we obtain 

1/( iV - 1 )/'c / 



(A.12) 



Therefore, the KN method developed for the nonlinear integral equation 
(A.l) is a special case of the GRS approach. 

We shall illustrate a numerical convergence of the iteration procedures 
(A.4)-(A.5) for an example [37] of the following nonlinear integral equation, 

V{x) = \ [\y(t)] 2 xtdt + l. (A13) 
2 Jo 

The exact solution of this equation has the form 

y(x) = l + ax, 

where a is given as solutions of a quadratic equation. There are two solutions 
of Eq. (A.9) with a « 0.40589 and a « 4.92744. 

For case of y^ = 1 + a^x, the iteration procedure (A4-A5) reduces to 

a m = [ i + i a (»-.) + i (a <»- 1 ) )2][1 __^ (1/3 + 2 l^»)] 

+ S nA ra « <N - 1, (| + ^). ( ' 
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Starting with = 0, let us calculate the differences defined as 

A (iV) = Uexact U (A16) 
O'exact 

where 

a exact — g — ~~ ^) ^ • 

The results are summarized in the Table I. 

From Table I, we can see a super-fast convergence of the procedure (A4- 
A5) for the nonlinear integral equation (A. 13). 

Starting with values of in range of < < 10, we have calculated 
the corresponding values of a^ N \ The results are summarized in the Table 
II. We note that in general the convergence of iterational procedure depends 
on initial step. If is chosen sufficiently close to the exact value we have 
a super-fast convergence. 
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Table I 



N 





1 


2 


3 


4 


AW 


1 


7 x io~ 2 


5 x io~ 4 


2 x 10- 8 


5 x 10" 17 



Table II 





aW 




a< 3 > 


a< 4 > 







0.37500 


0.40568 


0.40589 


0.40589 


0.40589 


1 


0.30000 


0.40352 


0.40588 


0.40589 


0.40589 


2 


-1.50000 


-0.03000 


0.37066 


0.40562 


0.40589 


3 


10.50000 


6.90957 


5.39040 


4.96679 


4.92778 


4 


5.25000 


4.94758 


4.92753 


4.92744 


4.92744 


5 


4.92857 


4.92744 


4.92744 


4.92744 


4.92744 


6 


5.10000 


4.93356 


4.92745 


4.92744 


4.92744 


7 


5.42308 


4.97200 


4.92787 


4.92744 


4.92744 


8 


5.81250 


5.05190 


4.93069 


4.92745 


4.92744 


9 


6.23684 


5.16756 


4.93897 


4.92747 


4.92744 


10 


6.68182 


5.31072 


4.95522 


4.92761 


4.92744 



15 



References 



[I] S. N. Bose, Z. Phys. 26, 178 (1924); A. Einstein, Sitz. 
Preuss Akad. Wiss. 1924, 3 (1924). 

[2] A. Griffin, D. Snoke, and S. Stringari, Bose-Einstein Condensation 
(Cambridge, New York, 1995). 

[3] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, 
and E. A. Cornell, Science 269, 198 (1995). 

[4] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, 
Phys. Rev. Lett. 75, 1687 (1995). In this work, the 
condensate could not be directly observed. See C. C. Bradley, C. 
A. Sackett, and R. G. Hulet, Phys. Rev. Lett. 78, 985 
(1997). 

[5] K. B. Davis, M. -O. Mewes, M. R. Andrew, N. J. van Druten, D. S. 
Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 
3969 (1995). 

[6] M. -O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, 
D. S. Durfee, and W. Ketterle, Phys. Rev. Lett 77, 416 
(1996). 

[7] D. S. Jin, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. 
Cornell, Phys. Rev. Lett. 77, 420 (1996). 

[8] M. -O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, D. S. 
Durfee, C. G. Townsend, and W. Ketterle, Phys. Rev. Lett. 77, 
988 (1996). 

[9] J. R. Ensher, D. S. Jin, M. R. Matthews, C. E. Wieman, and E. A. 
Cornell, Phys. Rev. Lett. 77, 4984 (1996). 

[10] C. J. Myatt, E. A. Burt, R. W. Ghrist, E. A. Cornell, and C. E. Wieman, 
Phys. Rev. Lett. 78, 586 (1997). 

[II] D. S. Jin, M. R. Matthews, J. R. Ensher, C. E. Wieman, E. A. Cornell, 
Phys. Rev. Lett 78, 764 (1997). 

[12] Mark Edwards and K. Burnett, Phys. Rev. 51, 1382 (1995). 

[13] P. A. Ruprecht, M. J. Holland, K. Burnett, and Mark Edwards, Phys. 
Rev. 51, 4704 (1995). 

[14] G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1966). 

[15] Mark Edwards, D. J. Dodd, C. W. Clark, P. A. Ruprecht, and K. 



16 



Burnett, Phys. Rev. A 53, R1950 (1966). 
[16] F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996). 
[17] A. Griffin, Phys. Rev. B 53, 9341 (1996). 
[18] A. L. Fetter, Phys. Rev. A 53, 4245 (1996). 

[19] R. J. Dodd, Mark Edwards, C. J. Williams, C. W. Clark, M. J. Holland, 
P. A. Ruprecht, and K. Burnett, Phys. Rev. A 54, 661 (1996). 

[20] F. Dalfovo, C. Minniti, S. Stringari, L. Pitaevskii, Phys. Lett. 
A227, 259 (1977). 

[21] Y. Kagan, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 
76, 2670 (1996). 

[22] Mark Edwards, P. A. Ruprecht, K. Burnett, R. J. Dodd, and C. W. 
Clark, Phys. Rev. Lett. 77, 1671 (1996). 

[23] E. M. Wright, D. F. Walls, and J. C. Garrison, Phys. Rev. Lett. 77, 
2158 (1966). 

[24] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[25] V. M. Perez-Garcia, H. Michinel, J. I. Cirec, M. Lewenstein, and P. 
Zoller, Phys. Rev. Lett. 77, 5320 (1996). 

[26] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 54, R4633 
(1996). 

[27] B. D. Esry, Phys. Rev. A 55, 1147 (1997). 

[28] A. L. Fetter, J. Low Temp. Phys. 106, 643 (1997). 

[29] V. L. Ginsburg and L. P. Pitaevskii, Sov. Phys. JETP 7, 858 (1958). 

[30] L. P. Pitaevskii, Sov. Phys. JETP 13, 451 (1961). 

[31] E. P. Gross, Nuovo Cimento 20, 454 (1961); J. Math. Phys. 4, 
195 (1963). 

[32] A. L. Fetter and T. D. Walecka, Quantum Theory of many-particle 
system (McGraw-Hill, 1971). 

[33] V. L. Ginzburg, Physics-Uspekhi 40, 407 (1997). 

[34] L. Vazques, L. Streit and V. M. Perez-Garcia, Nonlinear Klein-Gordon 
and Schrddinger Systems: Theory and Applications (World Scientific, 
(1996); M. J. Allowitz and P. A. Clarkson, Solitons, Nonlinear 
Evaluations and Inverse Scattering (Cambridge Univ. Press, New 
York, 1991). 



17 



[35] A. L. Zubarev, Sov. J. Part. Nucl. 4, 188 (1978). 

[36] E. Gerjuoy, A. R. P. Rau, and L. Spruch, Rev. Mod. Phys. 725 
(1983). 

[37] L. V. Kantorivich and G. P. Akilov, Functional analysis 
in Normed Spaces (MacMillan, New York, 1964), Chapter 18. 

[38] I. W. Burbridge, J. London. Math. Soc. 23, 135 (1948). 



18 



